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SUMARY OF ACTIVITIES 


During tha pe ri od of this grant, the authors succeeded in inpleraenting 
four out of the required six components necessary to define the spa t ial 
resolution requirements of candidate spa o eb o me micrcwave imaging systems. 

A brief discussion of these components follows; more complete discussions 
are contained in the conference publication in Appendix A and the draft 
technical publication in Appendix B. It is noted that precipitation 
measurements (e.g. , surface rain rates and water density and phase 
profiles) p lace the most severe demands of ary geophysical observable on 
the sp a tia l resolution, accuracy, and spectral coverage of a microwave 
imaging system [Stutzman and Brown, 1991] . Thus, we have targeted accurate 
and timely precipitation measurements as constituting the driving 
requirenents of such sensors. Our choice of maritime (as opposed to 
continental) precipitation is based on the greater utility of a spaoebome 
precipitation measurement system in observing the Earth over relatively 
inaccessible oceanic regions rather than land. 

( 1 ) qnwHai micrcphVBioai clcud and rain data. In order to derive 

realistic assessments of the perfo rm a n ce of c an d id a t e sensors, 
statistically meaningful cloud and raincell data has been obtained. Three 
sou rce s of data are used: (A) Simulated three-dimensional microphysical 

data from the Goddard Cumulus Ensemble model [Tao and Sinpson, 1889] . This 
consists of five frames sampled throu^icut the lifetime of a convective 
squall. The prevailing conditions are based on observed data from the 
flinhoi Atmospheric Research Program's Atlantic Tropical Experiment (GATE) . 
(B) A volume weather radar scan of an isolated multicell system observed 
during the Cooperative Huntsville Meteorological Experiment (OCHMEX, 1986) 
[Gasiewski, 1989]. (C) Synthetic cloud and raincell data based upon 
published spat ial and tenporal statistics of rainfall. 

(2) Forward radiative t ransfer model. Brightness maps for the above 
data at vir tual ly all microwave channels considered to be useful for 
tropospheric precipitation sensing have been ccnputed using the iterative 
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radiative transfer model described by Gasiewski and Staelin [1990]. these 
T constitute a set of "full resoluticn" mages fra. which particular 
chacnels subsets can be selected to study the performance of a candidate 
sensor, the channels under ocnsideraticn are at mcrovave wmta, 
frequencies fra. 6 to 410 GHz (9 channels), near the 22.235-, 183.31-, and 
325.15-®* water vapor lines (10 channels) , within the 5-- oxygen band (5 
channels) , an! near the 118-®z oxygen line (6 channels) . 

( 3 , ~in rattem ocnvolUtigBi spot patterns for ^itfractioh 

limited circular apertures are generated to preside the ooltispectral point 
spread functirh that would be expected fro. a candidate sensor. Each 
charnel is assumed to be diffraction limited, so that the miltispectral 
spot pattern exhibits successively decreasing 3-dB spot sizes as frequency 
increases. Currently, the degree of aperture illumination taper is 
selectable so that the effects of varying aperture efficiency and taper on 
the resulting retrieval accuracy can be studied. The candidate systems used 
in the siaulaticns are the three systems referred to in Appendix B. the 
BOS-B MMR I£0 system with 1.6- and 4.4-m apertures, and the large Spa 

Antenna (ISO niim-channel (6-410 GHz) and eight-charnel <«-«<“> 
systm. with IS- and 40-m apertures. The sutaillimet*: wave channels in the 
ni^chanrml ISA GEO system could potentially provide adeqiate spatial 
resoluticn for rainoell making using only moderately-sized apertures 

r Gasiewski, 1992]. 

The mlltispectral spot patterns for a candidate sensor are oonvolv 
with the full resolution imagery to obtain the antenna tenperature imagery. 
This imagery is subsequently saapled at the Nygiist resolution, effectively 
simulating the antenna scan process. After convolution, an ajpropria 
pseudorandom imtrvment error map is added to similate “ 

the detected signal. The noise level i- det^mhmd ^ ^ 
time, the receiver noise temperature, and the charnel bandwidth. Since 
these quantities can vary considerably among candidate dei gns, a 
representative set of receiver characteristics 

BOS KD* low-Earth orbiting system «d the ISRC large Space Antenna (ISA) 
geosynchronous system was chosen. 
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(4) Kartwnen-tegv? image SNR *"fllYBggi The Karhunen-Ioeve (KL) 

transform was us e d to estimate the number of observable spectral degrees of 
freedom in the noisy, convolved (received) imagery. Upon prewhitening the 
received imagery (so as to irrterccnpare channels on a noise-equivalent 
basis) , the KL transform was confuted. The eigenvalues of the KL transform 
are the signal -to-noise ratios (SNR's) of the unoarrelated spectral ™orUac 
of the imagery. Ccnpared to the full-resolution case, the MIMR imager y 
exhibits reductions of 1 to 3 dB in the SNR's of the second- and third- 
most dominant mode s for 4.4-m and 1. 6-m diameter apertures. Reductions in 
SNR for second- and third-most dominant modes in the ISA 40-m and 15-m 
systems are larger (up to 10 dB or more) . The dominant wnttoa in t hese 
systems exhibit no significant SNR reductions from the full resolution case 
to either of the respective smaller apertures. 

However, the interpretation of the KL modes and the SNR's 

for various spa t ial resolutions cannot be considered to be acnplete until 
these mode are incorporated into the precipitation parameter retrieval 
algorithms, discussed in ocnpcnent 6. 

(5) Mil tl snectral deoonvoiuHnn srndiea, Although not yet implemented, 

the in t a im a l 1 ate step of multispectral deconvolution is expected to 
simplify the resulting pre c i p ita t i o n retrievals by ( i-i rH p rj •the 

Bultispectral brightness image to the best d e g r e e possible. In the 
deconvolution, there will inevitably be a tradeoff between the resulting 
spatial resolution and noise in the deconvolved imagery. It is the 
manifestation of this tradeoff on area-averaged precipitation retrieval 
error which is of prime interest. The imagery from compone nt s 2-3 will be 
useful here in determining optimum deconvolution filters has**i on the 
sp at ial and spectral statistics of the underlying full-resolution imagery. 

(6) Area-averaged rain rate retrie val error nTYlIrrnfT. This final and 
most important ocnpcnent of the simulation will (upon complete 
Impl e m e n tation) produce surface rain rate, freezing altitude and water 
density and phase profile esti, mates hawed on the received, deconvolved 
imagery. 

TVro nonlinear statistical schemes, both optimal in their use of the 
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received data, are currently being implemented: (A) the single-step 

ronl ir»aar-«tatistical wtrie-ml eethcd, Uhl* uses KL rank reduction 
caduoe the mcbar of Up* channels prior to estimating eurfaoe rein rate 
art water density, art (B) the iterative statistical sethod, whld; uses the 
forward rrtiative tramfer relationships to determine successive linear 
statistical corrections to the retrieved parameters. Both of these settcds 
have been previously docunented [Gasiewski and Staelin, 1989; Kbo, 1 ]• 
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acMcnsicNS mo funs for future ncrk 


Oonclusicra from the first four of the six compo nents of the simulation 
study, including examples of full resolution, convolved, and Karhunen-Ineve 
imagery are contained primarily in the r epor t , "Simulations of the Effects 
of Spatial Resolution on Passive Microwave Remote Sensing of 
Precipitation' ' , in Appendix B. Although the quantitative effects of 
aperture size, main beam efficiency, channel selection, and antenna gain 
perturbations on simulated precipitation retrieved accuracy cannot yet be 
given 1 , several observations concerning effects on the received imagery and 
associated KL imagery can be made. 

An important feature is that the KL eigenvectors are relatively 
independent of storm type, nor do they depend critically on the stage of 
the raincell. In addition, the eigenvectors change only marginally as the 
aperture diame ter (i.e. , system resolution) changes. These observations 
suggest that the dominant inodes of spectral variation in microwave 
precipitation imagery are effectively unique far a given set of microwave 
charnels. However, as spatial resolution is lost, the eigenvectors 
increasingly esphasize the higher frequency channels. 

The convolved imagery and KL mode signal -to-noise (SNR) analysis 
unambiguously show a loss in spectral variance as aperture size is reduced. 
Far the apertures considered, the size reductions usually reduce the SNR in 
the dominant mrrto by only one or two dB, but reduce the higher order modes 
by much more (up to 10 dB) . The number of degrees of freedom in the 
p pfrrH p i m imagery varies far the different systems 
analyzed, but is approximately three to five. Although the effect that the 
KL vodm hsve on the accuracy of precipitation parameter retrieval is still 
to be the loss of information as aperture size is reduced is 

discernible. 

In microwave radiometry, antenna main beam efficie n ci e s of at 90-95% 
are considered to be r equir ed if deconvolution is not used. Note that this 


1 Future work on the problem of passive mul tiap ectral microwave 
imaging w ill continue beyond the end of this grant under NASA Graduate 
Student Research Program (GSRP) grant NCT-50903. 
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is the nall-to-null efficiency, as 3-dB efficiencies greater than 50-55% 
cannot practically be achieved. Such main beam efficiencies will require 
tapered aperture illuminations, as uniform tapers cannot yield main beam 
efficiencies greater than 86%. However, tapering necessarily reducing the 
system resolution. If the gain pattern is known precisely enough, seme 
deconvolution can be performed. By modelling the effects of antenna gain 
uncertainty as additional receiver noise (as shewn in Appendix B, section 
4) , the tradeoffs between the gain uncertainty, calibration uncertainty and 

illumination taper can be studied. 

Some oonments must be made about the spatial sanpling rate of microwave 
imagers. Wiile it is ccnmonly thought that sanpling in spatial increments 
to the 3-dB footprint of the antenna is adequate for passive 
microwave imaging this is, in fact, at least two to four times larger than 
the Nyquist sanpling increment. Thus, there is no guarantee against 
aliasing in any system that sanples at the «3-dB" rate. Of course, Nyquist 
sanpling requires significantly higher data rates and potentially higher 
nmcfcanical slew rates than 3-dB sanpling. One of the tradeoffs to be 
studied using the FWflRS will be the consequences of aliasing in sub-Nyquist 
sanpled imaging systems. 

Fu tur e simulation studies will utilize the KL modes in nonlinear 
precipitation parameter retrieval methods. The first method, currently 
being developed, will use a nonlinear single-step estimator which 
correlates the KL nodes with rain rate. Subsequently, an iterative scheme 
using a linear statistical operator at each step will be investigated as a 
means of inproving the performance of the single step nonlinear operator. 
©ie retrieval operators will incorporate spectral-dcmain nultispectral 
deconvolution operators (Rosenkranz, 1978] to improve the resolution of the 
brightness maps prior to inversion for the following parameters: surface 
rain rate, water density at five levels, mean ioe particle size, and 
altitude of nucleation. The analyses will attenpt to address 1.6-m and 4.4- 
m gbo systems as well as the three systems addressed in Appendix B. 
Additional three dimensional rainoell microphysical data will be included 

as available. 
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Ahatract-A numerical limulation of the latejlite-based mul- 
tiipectral ^uive microwave mapping and precipitation retrieval 
prorcs* is described. The simulation contains three major ^conv 
nofienti: (1) the forward radiative transfer calculation. (2) the 
ES observation simulation, and (3) the precipitation retrieval 
algorithm Secondary components include a karhunen-Loe e 
transformation module and antenna pattern deconvolution mo<b 
uI* tK simulation facilitates the investigation of ran I .late 
precipitation retrieval algorithm* using suIvNycpiist brightness 
femperature imagery and the comparison of the rela i tve me Ms 
of nrouosed radiometric svsteins such as the LOS Mullispcrlrai 
Imaging Microwave Radiometer (M1MR), the TRMM Microwave 
Imager (TRMM TM1). and future geosynchronous satellite sen- 
•ora. A procedure for retrieving rain rates using a statistical 
iterative approach is outlined. 

| Introduction 

In satellite-based passive remote sensi »"*' b ?^ 'll 1 '?"™ 
and infrared frequencies are commonly used to olrwne proper 
ties of atmospheric hydrometeors. Ail important fcatuic of in- 
frared imagery is the available high spatial resolution. However, 
diTu.ThV7arge extinction at infrared wavelengths such sen- 
JoT. Se unable'to probe through most clou, I cover. Relative to 
infrared aensing observations, mnltifrrquenry microwave 'ending 
provides a valuable coniplcinciitai rv capabi ity m being able ^to 
Srobe through clouds. This is nowiblc to * degrw winch depends 
on the particular frequency of observation and the hydronieteor 
density and sir.e distribution. For example, frequencies belo 
^ 6 Gill respond significantly to only very strong precipitation, 
while frequencies above - 220 Gilt respond to even light non- 
precipitating clouds such as cirrus. 

The relative simplicity of the forward transfer relationship 
at microwave frequencies along with the wide available tang of 
optical depths suggests that profiling of various precipitation pa- 
rameters (e.g. rain rate and hydrometeor density, siw, ami phase) 
m!ght*be possible. However, the non-linea, relation.l , .|» ^tween 
area- averaged precipitation parameters and the observed brig it- 
IS temperatures Is a major cause of error in retrieving these 
parameters. Impediments to improving retrieval accuracy in- 
clude: (1) sub-Nyquist spatial sampling caused by wide an enna 
beamwidths, (2) radiometric instrument noise, and (3) antenna 
U p.ll.rn nnccTl.inly. Since .!,, n.otogn, 
cosu for overcoming these impediments arc high, it is ot v ital 
importance to determine the available performance of candidate 
spaceborne sensor systems for cost /benefit tiadeoff pmpo. 

To aasess the relative merits of various P^rivT microvvave 
precipitation measurement systems, numerical simulations of h 
observation and retrieval process are being l^^^e * 
these simulations is to develop precipiUti°iM'* * l ' , ,< 1 

algorithms optimised for the channel, noise and gam pat n 
characteristics of individual sens™* lasting systems bniuj 1 st u< - 
led include the DMSP Special Sensor Microvvave/lmagei (SSN / ) 
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and the S|iecial Sensor Microwave/Temperalure Sounder (SSM /T- 
2); proposed sensors being investigated include the M1MK and 
the TRMM TM1. 

The simulation software, called Hie Passive Microwave Map- 
ping and Retrieval Simulator (PMMRS). provides a flexible end- 
o-ond siiiuilatiou of a mullisivectral passive microwave imaging 
and retrieval system. The PMMRS models three primary pro- 
cesses the forward radiative transfer process, the observation 
process, and the precipitation retrieval process. As shown in 
Fiem«- I. the forward radiative transfer process converts mirro- 
phvsical cloud data into upwelling brightness imagery, the ob- 
servation process transforms the upwelling brightness imagery 
into the brightness temperatures recorded by the antenna sen- 
sor while the retrieval process estimates atmospheric parame- 
ters (e.g., rain rale) and compares them to P*f* n J?‘"*"£ r,vod 
from the original microphysical cloud data. The PMMRS per- 
mits evaluating ret rievals for user-selected combinations of sensor 
antenna sped f^tiona. channel sets, noise levels, observation ge- 
ometry. and precipitation classes 

3 Forward R adiative Transfer Simula.liau 

Microplivsical cloud parameter data used in this investiga- 
tion is composed primarily of 3-dimensional synthesised data 
from the Cod, lard Cumulus Ensemble (GCE) simulation of a 
developing convective atorm sy stem [1). This data consists of a 
set of five time samples of a tropical squall devolo|>ed from ini- 
tial conditions measured during CAT E. \ertical profile* ° 
sure, temperature, relative humidity, li<|uid hydrometeor content 
and frozen hvdrometeor content are given for each pixel. Hor- 
izontal spatial resolution is 1.5 km. Vertical I resolution is typi- 
cally 1 km. Additional microphysical data based on a volume 
scan of a convective nnilticell system made by the C P-2 weather 
radar during the Cooperative Huntsville Meteorological Experi- 
ment (COIIMEX, 1986) is used. 

The microphysical cloud data wore mapped into upwelling 
brightness images (7' w ) using t lie iterative forward radiat ive trans- 
fer model developed by Gasiewski and Staclin [2]. I nis is a 
planar-st ratified scattering-baaed calculation that uses the spec- 
ified surface reflectivity and atmospheric parameters (trmpera- 
lure, pressure, relative humidity, liydromcleor distributions, etc. ) 
,t levels from Hie surface to ~ 25 km. Brightness temperature 
imagery over a calm ocean was computed for t.hr following mi- 
crovvave channels: 6. 10.60, 18.7. 23.8. 3G.5 50^3. >2.8. 53.5 )6. 
54.4. 51.01. 89, 118.75 ±(0.0, 0.12, 0.22. 0.3r. 0.C/. 1.2*, 2 0 '). 
166. 183.31 ±(0.0. 1.0, 3.0, 7.0), 220. 325.153 ±(0.0. 1.0. 3.0. , 0. 
9.0), 310. and 410 GHz. For example, Fig. 2a shows I lie lull- 
rcsolulion brightness temperature imagery for three oMhe six 
MINI It channels. The himodal structure seen m the 
the 18.7 Gllz channel it caused by cold ocean.c reflectance ul- 
sidi- the storm, a warming caused by thin absorbing cloud. 



the periphery of the storm. And a cooling due to ice scattering at 
the center of the storm system. 

3 Sensor Observation Simulation 

To simulate the observation process, four operations are per- 
formed: (I) an antenna gain pattern is calculated, (2) the up- 
welling brightness temperatures are convolved with the antenna 
gain, (3) the convolved imagery is sampled, and (4) pseudo- 
random instrument noise is added. A Bessel approximation is 
used in generating gain patterns for circular parabolic aperture 
antennas. The gain pattern generator requires the following user- 
supplied parameters: aperture diameter, illumination taper, fre- 
quency, and number of sidelobes to include in the gain pattern. 
The spatial response for each channel follows from the orbital 
parameters (e.g., satellite altitude). A two-dimensional convolu- 
tion of the calculated brightness imagery and the spatial response 
yields the antenna temperature imagery: 

T*(x,y) - ?(*,y)**Ta(x,y)*f n(x,y) (1) 

which is a blurred version of the underlying brightness imagery. 
The computed antenna temperature image for each channel is 
downsampled according to the instrument s highest effective spa- 
tial resolution. Gaussian pseudo- random noise is added to each 
pixel to simulate instrument observation error. The noise stan- 
dard deviation is Helermined for each channel based on I he avail- 
able bandwidth at that frequency, integration interval, and ex- 
pected svstem noise temperature. Figure 2b shows three channels 
of the MIMR imagery convolved with the sensitivity pattern of a 
circular 1.6-m linearly- tapered aperture at an orbital altitude of 
705 km. The 89-G)Iz channel has a 3-dB spot size of - 1.9 kin. 
and reveals most of the storm structure. At 6-GHz, the s|>ot size 
ls~28 km, and nearly all storm structure is lost . 


4. Karhuncn-Locve Analysis 

The original Tb and the convolved 7* imagery is analyzed for 
loss in observational degrees of freedom using a Karliuiien«lxx*vc 
(KL) transformation (also called principal component analysis). 
The unitary KL transform is used to rank-order the dominant 
spectra] modes in the 7* and Ta imagery by the mode signal- to- 
noise ratio (SNR) (3). The modes are determined by the eigenvec- 
tors of the spectral covariance matrix. The eigenvectors provide 
a measure of the contribution from each channel to the trans- 
formed image. Associated with each eigenvector is an eigenvalue 
which is the variance for the particular 1\L mode. To account for 
observational noise variations among the channels, prewhitening 
the brightness imagery prior to compulation of the covariance 
matrix is performed. 

The two most dominant KL images (i.c., those with the two 
largest SNR’s) for the six channel full- resolution MIMH imagrrv 
are shown in Fig. 3. The adjoining eigenvector plots reveal that 
the most dominant mode is essentially the spectral difference l>c- 
tweeii the 23.8- and 89-GHz channels. We interpret this as the 
strong scattering signature from ice aloft. The second KL mode is 
a combination of the scattering signature from the high frequency 
channels and the spectra) difference between the 6- and 10.69- 
GHs channels, and appears to be sensitive to the radiomctricallv- 
thin periphery of the storm. The KL-mode SNR’s for the full- 
resolution an 3 convolved MIMR imagery (Fig. 4) can be used to 
determine the loss in number of observable degrees of freedom 
caused by the antenna pattern’s broadness. If we arbitrarily as- 
sume that a KL mode contains useful information only if it has 
a SNR>2.0»3dB, the number of observable degrees of freedom 
decreases from approximately three to two u|K»i convolution of 
tile MIMR antenna pattern. 


5. Precipitation Retrieval Simulation 

This component of the PM MRS provides an environment for 
developing and optimising precipitation retrieval algorithms for 
multispectral sub-Nyquist sampled microwave imagery, lo take 
advantage of the spatial resolution available from high-frequency 
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channels, multispectral deconvolution of the antenna gain pattern 
is first performed. Two deconvolution modules are being devel- 
oped. Both methods use information from the highest resolution 
(i.e.. highest frequency) channels to improve I ho resolution of the 
lower frequency channels. The first method uses the well-known 
Linear Minimum Mean Squared Error (LMMSE) approach in the 
spatial domain, such that: 

t fl = T5T a (2) 

where T * is a vector consisting of all pixels from all channels 
having any significant correlations with the brightness at the lo- 
cation of interest. Such correlations can be due to either rain cell 
structure, spectral scattering or absorption structure, or antenna 

pattern filtering. The determination matrix 25 is given by: 

T5 - T^t 9 t a (Rt a t a + TLnV' ( 3 ) 

Tiij is the joint spatial-spectra! covariance matrix between the 
multispectral 2-dimcnsional signals t and j. 

The second deconvolution method, developed by Rosenkranz 
|4], uses a determination matrix 75(u,t>) in the spatial Fourier 
domain, ami thus assumes stationarity in the statistics of (he 

underlying brightness imagery. In this method Tb is found from: 

a(x, v)) = ^(u,v).F{T,t(T, 5 ()} (4) 

w here T is the 2-dimcnsional Fourier transform and u, r are spa- 
tial frequencies. The expression for 23(u,u) can be found in |4). 
To obtain T«. the inverse Fourier transform is used. 

Subsequent to deconvolution, rain rate retrieval algorithms 
are to be implemented. Parameter estimation techniques under 
consideration include LMMSE algorithms, Wiener filtering, and 
non -linear statistical and iterative estimators. For optimum re- 
trievals, the non-linear relationship between brightness temper- 
atures and cloud and precipitation parameters must Ik* consid- 
ered. Here, a statistical iterative retrieval technique baaed on 
the piecewise- 1 inear relationship between incremental changes in 
brightness temperature and incremental changes in surface pre- 
cipitation rate (via the absorption and scattering coefficients) is 
of interest. This technique is similar to the LMMSE technique 
except that 7J becomes an incremental predictor matrix. and 
A7 'b and A/f (where AT? is the incremental change in rain rate) 

replace Ta and Tb, respectively, in Eqs. 2. 3, and 4. The sta- 
tistical iterative technique has been proved useful in retrieving 
relative humidity profiles [/»]. 

6. Piacumu 

At this time both the forward radiative transfer and the 
sensor observation comj>onents of the PM MRS have been imple- 
mented. In addition, a study of the number of observable degrees 
of freedom for sensors systems with different antenna aperture di- 
ameters and observation geometries has been performed. Some 
results from this study arc in Table I. Although we cannot yet 
quantify the available retrieval accuracy of candidate systems, the 
loss in number of observable degrees of freedom as aperture size 
is reduced is readily seen. Future implementation of the deconvo- 
lution and parameter retrieval components of the PM MRS will 
improve our ability to determine (quantitatively) the available 
accuracy of existing and future spaccbor no passive microwave 
meteorological systems and lo design optimized sciisoi-*|>ccific 
multispectral retrieval algorithms. 
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Figure 3: The two most dominant Karhunen-Loeve modes for 
the full-resolution imagery in Fig. 2a, along with the associated 
eigenvectors. 




Figure 2: Representative MIMR brightness temperature imagery 
for G, 18.7, «md 89-Gllz: (a) full resolution. and W »fier convo- 
lution with the gnin pattern of ft circular 1.6-m linearly-tapered 
aperture. The respective 3-dB apot lizet are 28, 9 and 1.9 km. 
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Figure 4: SNR ranking for full- resolution and convolved MIMR 
imagery. 


Table I: Number of Observable Degrees of Freedom 


Sensor 

Aperture 
diameter (m) 

# 

channels 

# degrees 
of freedom 


00 

9 

7 


30 “j 

9 

4 

mTmk 

oo 

G 

3 


1.6 

G 

2 


Notes: 

(1) Assumes circular liiicarly-taporcd aperture field. 

(2) D * oo implies full resolution imagery. 

(3) Degrees of freedom observable only if SNR>2.0=3dB. 
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Ab«aa-A numerical .io.ul.lion of the satellite-b»s«l mult, spectral passive microwave 
precipitation mapping and retrieval problem i. described. The purpcwe of .hi. simulation ., .o 
provide uttenn. and radiometer .y.tem. engineer, with ,»antita.i.« performance pred,c..on. 
for proposed .paceborne passive microwave imaging tensors. The sunulat.on contains s« ma- 
jor components. (1) microphysical cloud and raincell data compilation, (2) forward radiative 
transfer calculations, (3) sensor observation simulations, (4) Karhunen-Loeve transformation 
and degree of freedom calculations, (5) antenna pattern deconvolutions, and (6) prec.pitat.on 
parameter (e.g., rain rate) retrieval algorithm development. The simulations facilitate .nvest, ga- 
llon of precipitation parameter retrieval using low-pass filtered brightness temperature imagery 
from candidate passive microwave systems. The investigation also facilitates comparison of the 
relative merit, of propo«d radiometric systems such as the EOS Multispectral Imaging M.- 
crowave Radiometer (MIMR) and geosynchronous satellite sensors. Initial simulation resuHs for 
components # 1-4 are described herein, and procedures for deconvolution (component # 5) and 
a statistical iterative approach to rain rate retrieval (component # 6) are outlined. 
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1. Background 


In satellite-based passive remote sensing, both microwave and infrared frequencies are 
used to observe properties of atmospheric hydrometeors. An important feature of infrared im- 
agery is the high spatial resolution available from instruments of moderate aperture size (< 1 m), 
even at geosynchronous distances. Due to the large hydrometeor extinction at infrared wave- 
lengths, such sensors are unable to probe through most cloud cover. In contrast, multifrequency 
microwave remote sensing provides a valuable complementary capability in being able to probe 
through clouds and even some precipitation while retaining useful sensitivity to hydromete- 
ors [1, 2]. The degree of sensitivity depends on the particular frequency of observation, the 
hydrometeor phase (liquid or ice) and the hydrometeor density and size distribution [3]. For ex- 
ample, frequencies below ~ 6 GHz respond significantly to only very strong precipitation, while 
frequencies above 220 GHz respond to even light non-precipitating clouds such as cirrus. 

Compared to the infrared, the precipitation probing capability of microwave observa- 
tions greatly facilitates the measurement of precipitation parameters such as surface run rate. 
Precipitation sensing capabilities have been demonstrated by a number of airborne and space- 
borne microwave instruments, for example, the Advanced Microwave Precipitation Radiometer 
(AMPR) [4], the Millimeter-wave Temperature Sounder (MTS) [5], the Nimbus-E Microwave 
Spectrometer (NEMS) [6], the Scanning Multichannel Microwave Radiometer (SMMR) [7], and 
the DMSP Special Sensor Microwave/Imager (SSM/I) [8]. Measurements of surface precipita- 
tion rate facilitate hydrological studies involving water budget and runoff analyses.ln addition, 
the relative simplicity of the forward radiative transfer relationship at microwave frequencies 
along with the wide available range of cloud penetration depths suggest that profiling of some 
precipitation parameters (e.g. rain rate and hydrometeor density, size, and phase) might be pos- 
sible [9]. Measurements of water density profiles in precipitating systems would be particularly 
useful for severe storm tracking and estimating atmospheric heating profiles, which are in turn 
important in understanding global heat transport [10]. 

For applications in hydrology, severe storm tracking, operational meteorology and global 
change study, area-averaged estimates of the forementioned precipitation parameters with spatial 
resolutions of ~ 1-5 km and observed at time intervals of ~ 30—60 min are generally considered 
to be ultimately desirable [11, 12]. Further improvements in spatial or temporal resolution would 
yield only small marginal benefits. Degradations in spatial resolution to — 25 - 50 km would still 
yield useful data for some of the above applications, although the utility falls off quickly beyond 
this range. However, the accuracy of area-averaged estimates depends strongly on the spatial 
resolution, radiometric sensitivity and channel set of the observation system. Indeed, the forward 
transfer relationship is decidedly non-linear with respect to most hydrometeors parameters, 
especially at frequencies above ~ 23 GHz. In addition, there is often no unambiguous inverse 
relationship, particularly if observations are made over only a narrow range of frequencies. These 


14 



two effects, combined with (1) low-pass spatial sampling caused by large antenna footprints 
relative to the spatial scales of precipitation, (2) radiometric instrument noise, and (3) antenna 
calibration and gain pattern uncertainty compromise the accuracies obtainable from passive 
microwave imagery. Since the costs for overcoming the technological impediments causing these 
problems are quite high, it is of great importance to be able to predict the performance of 
candidate spaceborne sensor systems prior to implementation for cost /benefit tradeoff purposes. 

During the last decade, several new geosynchronous (GEO) and low-Earth orbitting 
(LEO) passive microwave precipitation sensors with enhanced spatial resolution capabilities 
have been either studied, proposed or recently deployed: 

1. Microwave Atmospheric Sounding Radiometer (MASR) GEO system: ~ 4.4 aperture [13] 

2. Geosynchronous Microwave Precipitation Radiometer (GMPR): ~ 4 aperture [14] 

3 ] LaRC Large Space Antenna (LSA): ~ 15 - 40 meter GEO imaging sensor [12] 

4. DMSP Special Sensor Microwave/Temperature Sounder (SSM/T-2) 

5 ] Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) [15] 

6. Multispectral Imaging Microwave Radiometer (MIMR): ~ 1.6-m LEO system for deploy- 
ment during EOS-B [16] . 

7. Jet Propulsion Laboratory Synthetic Aperture Imaging Radiometer (SAIR) [17] 

This list is not exhaustive, but simply representative. For the first three (GEO) systems, the 
most important advantage stems from the high temporal resolution available from the geosyn- 
chronous vantage point, although large apertures and/or high frequency channels are required to 
obtain reasonable spatial resolution. For the other systems, the low Earth orbitting configura- 
tion considerably simplifies the structural problems (and expenses) associated with large space 
antennas. However, the advantage of high temporal resolution is lost due to long revisit times. 
It can only be regained by implementing a small fleet of LEO sensors. 

To assess the relative merits of the above listed and other candidate passive microwave 
systems for precipitation measurement, numerical simulations of the multispectral imag- 
ing and retrieval processes as applied to spaceborne passive microwave precipitation measure- 
ment systems are being performed. The simulations are being based on realistic three-dimensional 
geophysical models, accurate sensor models, and optimal precipitation parameter retrieval algo- 
rithms. This paper discusses recent progress and the current state of the numerical simulations. 
A goal of these simulations is to develop precipitation parameter retrieval algorithms optimized 
for the channel set, noise and gain pattern characteristics, and observation geometry of a par- 
ticular sensor. 

The overall scope of the simulations includes all of the proposed sensors listed above, 
although the immediate application is being directed to three specific systems: (1) EOS MIMR 
(#8), using six channels at 6, 10.69, 18.7, 23.8, 36.5 and 89 GHz from low Earth orbit (705 km 
altitude), (2) LaRC LSA (#3), using nine window channels at 6, 10.69, 18.7, 36.5, 89, 166, 
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220, 340 and 410 GHz from geosynchronous orbit, and (3) LaRC LSA (#3), using two window 
channels at 18.7 and 36.5, one water vapor channel at 23.8 GHz, and five oxygen channels at 
50.3, 52.8, 53.596, 54.400 and 54.940 GHz (a total of eight channels) from geosynchronous orbit. 

Although we have targeted precipitation as the geophysical process of greatest interest, 
the retrieval of several other geophysical processes might be considered. For example, spaceborne 
measurements of temperature and water vapor profiles are also viable using passive microwave 
systems. However, the spatial and/or temporal resolution needs for measurement of these vari- 
ables are not as demanding as for precipitation or cloud water. When compared with water 
vapor or temperature, precipitation parameter retrievals are of greater concern for two primary 
reasons: (1) the relationships between the parameters and the observables are highly nonlinear 
and even bimodal, and (2) convective precipitation cells and other storm features (e.g., hurricane 
eyewalls) can be as small as a few kilometers in size, and typically a factor of five smaller than 
the spatial scale of most temperature or water vapor structure. Thus, if practical precipitation 
measurement systems can be identified, temperature and water vapor measurements using the 
same systems will not require significantly more sophisticated hardware or retrieval algorithms. 

The same can be said for sea ice detection. Here, spatial resolutions smaller than one 
kilometer are desirable, but the detection mechanism is relatively straightforward compared to 
precipitation retrieval algorithms and the time evolution of ice movement is relatively slow. 
Essentially, the problem requires a set of channels that exhibit a sensitivity difference to ocean 
water and sea ice, but can still probe through the atmosphere. This must be done at least 
once a day (or so). Here, high frequency window channels (e.g., at 37, 90, 166 and 220 GHz, 
including both vertical and horizontal polarizations) can be used since sea ice occurs only in 
cold latitudes where water vapor screening is small and high cloud opacity is not persistent. 
The necessary spatial resolution is obtainable by virtue of the narrow antenna beams available 
at these frequencies using apertures of practical size. In contrast, precipitation measurements 
are best performed using low frequency channels (e.g., 6, 10, 18 and 37 GHz) for which sub- 
Nyquist res ol ution is the overriding technological impediment (and hence, the overriding cost 
issue). Moreover, a geosynchronous system is not practical for sea ice observation due to the 
large incidence angle near the polar regions. 

Soil moisture can evolve quickly, particularly during precipitation events. For hydro- 
logical purposes, spatial scales of several kilometers are ultimately desirable. However, since 
the soil moisture retrieval problem is predominantly linear [19), sub-Nyquist spatial resolution 
should cause no significant degradation in the resulting area-averaged soil moisture retrieval er- 
ror. Moreover, since frequencies for soil moisture measurement are relatively low (near L-band), 
spaceborne systems with adequate spatial resolution will likely be designed as single-channel 
thinned-aperture interferometers [18). Again, the spatial resolution problem as applied to mul- 
tichannel precipitation sensing is of greater immediate interest. 

t6 



In contrast to water vapor, temperature, sea ice andwil moisture, the radiative transfer 
through precipitation results in strong nonlinearities between the underlying Mrometeor pa- 
rameters (e.g„ water density, phase and particle site distribution) and the upwellmg brightness 
temperatures. For sub-Nyquist systems, these nonlinearities can cause ambiguities in area- 
averaged retrievals of precipitation parameters. For example, using 37-GHt data from - 
belt el el. [20], a 50-mm/hr rain cell occupying 10 percent of the area of a sensor s footprint 
corns approximately the same brightness as a 1-mm/hr brittle occurring over the entire area 
of the footprint, even though the arc-averaged rain rates differ by a factor of five. 

It is for the above reasons that we have focussed our simulation efforts on precipitation 
parameter measurements, with particular emphasis on maritime precipitation. As a practical 
matter, maritime precipitation is observable only from space, and moreover occurs over a larger 
portion of the Earth’s surface than continental precipitation. In addition, continental precipn 
tation is more readily amenable to ground-based radar observation, and ,s currently observed 
over a significant fraction of the industrialised northern hemisphere. Thus, we have given initial 
emphasis to simulations over ocean backgrounds. We have also emphasised geosynchronous sys- 
tems. Precipitation events (particularly the convective type) evolve swiftly, with typical lifetimes 
of one to three hours. Thus, it is of particular interest to simulate precipitation s 

from geosynchronous orbit where short revisit times (less than one hour) can be achieved [2 ]. 

The simulation software, called the Passive Microwave Mapping and Retrieval Simula- 
tor (PMMRS), provides a Hexible end-to-end simulation of a multispectral passive microwave 
imaging mid retrieval system. The PMMRS models three primary processes: the forward ra- 
diative transfer process, the observation process, and the precipitation retrieval process. As 
shown in Figure 1, the forward radiative transfer process transforms three-dimensional micro- 
physical cloud data into upwelling brightness imagery, the observation process transforms the 
upwelling brightness imagery into the brightness temperatures recorded by the antenna sensor, 
while the retrieval process estimates atmospheric parameters (e.g., rain * ““ 

to parameter, derived from the original microphysical cloud data. The PMMRS permits eval- 
eating retrievals for user-selected combinations of sensor antenna specifications, channel se s, 
noise levels, observation geometry, and precipitation types. The overall goal of these mapping 
simulations is to quantitatively assess the effects of multispectral blurringonarea-averagede, 
rors in retrieved precipitation parameters such as surface rain rate and liquid and ,ce density 
profiles In studying the resolution problem, the authors recognise that some severe mechanical 
Lid electrical requirements might be imposed in the implementation of the tmndidate imaging 
systems. However, the primary purpcwe of this study is to define the imaging needs ba«d on the 
gMphyrid phenomena of interest, and to secondarily consider the (albeit formidable) hardware 

issues. 
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2. Forward Ra diative Transfer Simulation 


Microphysical raincell parameter data used to date is composed primarily of 
dimensional synthesized data from the Goddard Cumulus Ensemble (GCE) simulation «t a 
developing convective storm system (22). This data consists of a set of five time samples tit a 
tropical squall computed using a numerical cloud evolution model from initial conditions *K|i- 
sored during the Global Atmospheric Research Program Tropical Atlantic Experiment (GAft|. 
Vertical profiles of pressure, temperature, relative humidity, liquid hydrometeor content a#fl 
frozen hydrometeor content are given for each of 64 x 64 pixels. The horizontal spatial r<§r 
olution is 1.5 km. The vertical resolution varies, but is typically 1 km, with twenty discrete 

levels. 

Additional raincell microphysical has been obtained from a volume scan of a convective 
multicell system made using the CP-2 weather radar during the Cooperative Huntsville Meteo- 
rological Experiment (COHMEX, 1986) [23). The reflectivities were converted to hydrometeor 
densities and mean size parameters using the Marshall-Palmer (MP) and Sekhon-Srivastava (SS) 
size relationships. Particle phase was assumed to be liquid below the freezing level, solid above 
the temperature of ice nudeation (assumed to be -30° C ), and linearly mixed in between. A 
0.5 dB ice reflectivity correction was applied [24]. 

The microphysical raincell data was mapped into upwelling brightness images (7 b) using 
the iterative forward radiative transfer model developed by Gasiewski and Staelin [25]. This 
is a planar-stratified scattering- based model that uses the spedfied surface reflectivity and at- 
mospheric parameters (temperature, pressure, relative humidity and hydrometeor densities, size 
distributions and phases) at levels from the surface to ~ 20 km altitude. Brightness tempera- 
ture imagery over a calm ocean was computed for the following microwave frequencies: 6, 10.69, 
18.7, 23.8, 36.5, 50.3, 52.8, 53.596, 54.4, 54.94, 89, 118.75+(0.12, 0.22, 0.37, 0.67, 1.27, 2.07), 
106, 183.31-K0.0, 1.0, 3.0, 7.0), 220, 325.153+(0.0, 1.0, 3.0, 7.0, 9.0), 340, and 410 GHz. These 
frequencies indude all significant microwave “window-channels” along with those channels near 
the water vapor and oxygen absorption lines that are deemed to be the most essential for tro- 
pospheric measurement purposes. For the absorption line channels, the offset frequencies from 
the respective line centers are listed in parentheses. The resulting sets of brightness maps are 
catalogued in Table 1. 

Typical full-resolution brightness imagery for the microwave window channels is shown 
in Figs. 2a-d. These figures include a representative set of three GCE time samples and the 
single COHMEX sample. Here, whiter areas correspond to warmer brightnesses. In all imagery, 
the warming caused by precipitation over the reflective ocean surface is readily seen in many of 
the channels below ~ 23 GHz. This is expected since precipitation is predominantly absorptive 
at these low frequencies. Thus, emission from the relatively warm hydrometeors increases the 
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upwelling brightness above the cold values characteristic of clear-air regions over ocean surfaces. 
In contrast, scattering of the cold cosmic background radiation by ice over the storm cores 
is apparent in all higher-frequency window channels (~ 18 GHz and above). The bimodal 
brightness structure clearly seen in the imagery of the 18.7- and 36.5-GHz channel is caused by 
cold oceanic reflectance outside of the storm, a warming caused by thin absorbing clouds near 
the periphery of the storm, and a cooling due to ice scattering at the center of the storm. 

Corresponding full-resolution imagery for the oxygen and water vapor channels for GCE 3 
is shown in Figs. 3 and 4. As expected, the storm structure is least visible in the most opaque of 
these channels, that is, those channels with greatest integrated opacity. In the case of the oxygen 
channels, the altitude of the raincell top can be estimated by virtue of the multiplicity of channel 
probing depths obtainable using a series of channels with successively increasing opacities (23). 
These high-resolution brightness maps and others for similar precipitation cases are used as the 
basis for all ensuing statistical calculations and imaging and retrieval studies. 
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Candidate sensors are modelled by aperture size, aperture illumination field distribution, 
orbital altitude h, and field-of-view (FOV). In addition, each sensor channel n = 1 ... TV is mod- 
eled by the following characteristics: center frequency /„, bandwidth W ny and double-sideband 
receiver noise temperature JYec«- A diffraction-limited aperture is assumed for all channels, 
thus giving rise to spot resolutions that increase with channel center frequency. Collectively, 
these parameters determine the single-spot brightness sensitivity of each channel: 



where 7sysn = 7recn + TanU 16 system noise temperature (receiver plus antenna), and r is 
the integration (or dwell) time for a single sample. For Earth observation, typical antenna tem- 
peratures r^tn between 100 and 300 K, depending on the center frequency and geophysical 
state. 


For GEO systems, r is derived from the image FOV, orbital altitude, and the total scene 
acquisition time T„ assuming uniform rectilinear sampling of the subsatellite (equatorial) scene 
at intervals determined by the Nyquist sampling criterion for the highest frequency channel 
(Fig. 5a). Denoting this sampling angle by Bs (in radians), we have: 


r -T 1***1’., 

at0 ~ 1 * "fov" 


where h = 3.5885 x 10 4 km and the FOV is in km 3 . The parameter r), is the scan efficiency, 
which is the fraction of time spent acquiring radiometric scene data as opposed to performing 
supporting operations such as calibration and scan turnaround. For the GEO simulation studies, 
the FOV, total acquisition time and scan efficiency are taken to be 1,000 km x 1,000 km, 
T, * 30 minutes and 85%, respectively. 


For LEO systems, r is derived from the angular 6wath width Bw and altitude h, assuming 
a cross-track raster scan with the subsatellite raster spacing determined by Nyquist criterion for 
the highest frequency channel (Fig. 5b): 

B%h I h + 6356.8 _ 

T,eo B w V 3 986 * 105 Vt 3 

where Bs, Bw are in radians and h is in km. For purposes of antenna characterization, we use 
a cross-track scanning geometry of angular width 90°, r) $ =45% and h = 705 km for the MIMR 
simulations. 1 We also assume that scene mapping takes place using single (non-redundant) 

1 While in view of the f&ct that MIMR is a conical scanning instrument, an intracomparison of antenna char- 
acteristics can be most easily accomplished using the chosen raster scan geometry. Conical scans would produce 
similar aensitivitiea and resolutions, hence a meaningful comparison can be made. 
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receivers for each channel. If redundant receivers are used to provide parallel scene mapping, 
the effective integration time can be increased by the number of redundant receivers AW- This 
reduces each ass n value by 1 / V A f rec • 2 

To simulate the observation process, four operations are performed: (1) an antenna gam 
pattern (or, spot pattern) is calculated, (2) the upwelling brightness temperatures are convolved 
with the corresponding spot pattern, (3) the convolved imagery is sampled at the Nyquist rate, 
and (4) pseudo-random instrument noise is added. Discretized spot patterns are generated from 
the antenna gain pattern sampled at appropriate locations on the Earth’s surface. Gam patterns 
are based on the specific aperture field distribution and wavelength. Without significant loss of 
generality, apertures are assumed to be circular of diameter D, and have a linearly polarized 
aperture field with uniform phase taper and radial amplitude taper: 

r(r,») = £.x|l-(|) 1 ]’’ (4) 

Here, p is a parameter that describes the degree of illumination taper. The resulting family of 
gain patterns are described by Bessel functions [26]: 


Sn(p;0) = 9o 


Jp+i (^»infl) 




P+1 


( 5 ) 


where g 0 is a constant, and c is the speed of light. By sampling the continuous gain function, 
the diagonal elements of the gain matrix g are obtained: 


h = 


ffi(p;0«j) 0 

o p2(p;^) 


9N{P\0ij) J 


( 6 ) 


( 7 ) 


On = tan * 1 (tV‘ j + jJ ) 

where A R is the horizontal spatial resolution of the brightness maps (see Table 1). 

Since the Fourier transform of the antenna gain is related to the aperture electric field 
correlation function [27], the Nyquist sampling interval can be determined from the aperture 
size. This follows from the fact that the aperture field correlation function is zero for spatial 
displacements greater than D. Accordingly, if no aliasing is to occur, the maximum angular 
sampling interval 6s must be: ^ 

( 8 ) 


6s = 


2 Df N 


1 Redundant receivers may be the only practical mean, of meeting the .lew requirement, for «ome large filled- 

aperture .y.tem., for example, the MIMR 4.4-m system. 
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The corresponding spatial sampling interval is: 


Rs - Ms 


ch 

2 T7d 


( 9 ) 


This is the Nyquist sampling interval. It must be distinguished from the the 3-dB spot size (or, 
“3-dB resolution”) of the highest frequency channel, even though the 3-dB measure is commonly 
referred to in discussing system resolution. Depending on the aperture illumination taper, Rs 
is from two to nearly four times smaller than the 3-dB spot size at the highest frequency, which 
follows: ^ 

^3dB,N = Mmb,N = = 2q (p)^s ( 10 ) 

where the constant a (Table 2) is a weak function of p, ranging between approximately one to 
two ([26], pp.195). Note that sampling at Rn does not guarantee that all spatial frequencies in 
the underlying brightness imagery are measured. 


The discretized spot pattern is subsequently truncated so that the convolution in Eq. 14 
is over a finite range of indices. To determine the indices at which the pattern is truncated, the 
beam efficiency functions i?b(p; 0) are used. These functions are the fraction of beam energy in 
the conical solid angle from the beam axis out to a given angle $ (Fig. 6): 



The curves show that 3-dB beam efficiencies greater than ~ 55% cannot be obtained using 
simple aperture tapers. Note that the illumination taper p affects the main beam efficiency (or, 
null-to-null efficiency) considerably, particularly for weak tapers (p = 0 or 1). Specifically, main 
beam efficiencies greater than ~ 86% cannot be achieved with a uniform taper (p = 0), but 
require a more gradual taper (p > 1). However, increasing p also increases the 3-dB spot size 
by virtue of the factor o in Eq. 9. The aperture efficiency: 

“ (F7T7 (12) 

is a measure of the effective use of aperture area in increasing main beam gain, and is also 
reduced (Table 2) as p increases. 

Using Eq. 11, we find that truncation near the third null for the lowest frequency channel 
insures that virtually all (> 99.8%) of the beam energy is accounted for in the discretized spot 
patterns for all channels and for all tapers p > 1. An minor exception occurs for the uniform 
taper (p = 0) for which only 96% of the beam energy is contained within the third null for 
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.he lowest frequency channel. To account for both (I) radiation from far sidelobes and (2) the 
undetermined constant s„ the truncated gain matrix is subsequently normalised so that: 


i»ig 

E 


_ 

E h = * 


•*-'g J*"7g 


(13) 


where J,, is rero for |.| > ig and |;| > ;g, and 7 is the identity matrix. This insures that the 
spatial response is unitary over the (2tg + 1) (2jg + 1) values. 

To simulate the observation process, discretized gain patterns are subsequently convolved 
with the brightness temperature imagery then corrupted by additive white noise. Mathemati- 
cally: 


i'*-oo ;'»-00 


(14) 


where: Tbh is the multispectral scene brightness at pixel i,j (i = l...»max,J - 1 •••>max), 
7a ■■ is the associated multispectral antenna temperature at pixel i,j, is the sensor s mul 
ti spectral discretized spot pattern (or, gain) matrix, and W,, is white Gaussian pseudo-random 
observation noise with covariance matrix: 




<4i 


0 

6 <rf 2 


Mj 


(15) 


T h 


where * is the transpose operator and 6 is the Kronecker delta function. In Eq. 14, - denotes 
two-dimensional convolution. In performing the convolution, the brightness images are padded 
with dear-air spectra at locations outside the range of available data, effectively extending the 
images beyond imaxJmax pixels. Note that since the map resolution A R is smaller than the 
resolution of the sensor’* sampling grid R s , the map noise standard deviation *r» must be 
adjusted relative to the single-spot standard deviation assn by the relative sampling rates. 


For the LSA eight-channel system, the map resolution A R is considerably smaller than the 
Nyquist sampling interval, although A R is comparable to R s for the MIMR and LSA nine- 
channel systems. Overall, the above simulation of the observation process is equivalent to 
sampling the raster at a resolution f?s, then resampling to a standard resolution A R. For e 
MIMR and LSA systems, the assumed system sensitivity parameters, resolutions and associated 
integration times and single-spot sensitivities are listed in Tables 3 and 4. 
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Figs. 7a-c show the MIMR six-channel imagery before and after convolution with the 
sensitivity patterns of circular 4.4-m and 1.6-m quadratically-tapered (p = 1) aperture distri- 
butions at a 705-km orbital altitude. Here, R s = 0.27 and 0.74 km (respectively) while the 
3-dB resolutions are R 3< j B jj = 0.68 and 1.89 km. For illustration, the 210- minute GCE frame 
(GCE 4) has been chosen, and instrument noise has not been added. In practice, the convolved 
antenna temperature imagery would actually be downsampled at the Nyquist resolution of the 
highest frequency (89-GHz) channel. To obviate the need for subsequent upsampling to the 
resolution AR of the unconvolved brightness maps, we show the antenna temperature imagery 
at this same resolution. 

With a 4.4-m aperture, the fine storm structure observed in the full resolution imagery is 
revealed at frequencies as low as 10.69 GHz, and the 6-GHz channel shows some of the overall 
shape characteristics. However, with the 1.6-m aperture (the size proposed for MIMR during 
EOS-B [16]), the fine storm structure is visible only at frequencies equal to or higher than 
18.7 GHz. At 6-GHz, the spot size for the 1.6-m aperture is ~ 28 km, and even the gross shape 
seen in the full resolution imagery is lost. In both convolved cases the 89-GHz channel has a 
small enough 3-dB spot size (~ 0.7 and 1.9 km for the 4.4- and 1.6-m apertures, respectively) 
to reveal virtually all of the features seen in the full resolution imagery. 

Figs. 8a-c show the nine-channel GEO LSA imagery for GCE 4 both before and after con- 
volution with the sensitivity patterns associated with circular 40-m and 15-m diameter apertures 
with quadratically-tapered (p = 1) field distributions. Even with the 40-m aperture, much of the 
fine storm structure is lost in the process of imaging. However, some of the bimodal brightness 
structure nesur the storm edge remains observable at 36.5 GHz. With a 15-m aperture, very little 
bimodal structure remains at even 36.5 GHz. Although the blurring at frequencies of 166 GHz 
and higher is small for both 15- and 40-m apertures (spot sizes are less than ~ 6 km), there is 
considerably less information on fine storm structure at these frequencies. Similarly, Figs. 9a-c 
show the eight-channel (low-frequency) GEO LSA imagery for GCE 4 both before and after con- 
volution with the sensitivity patterns associated with circular 40-m and 15-m diameter apertures 
with quadratically-tapered (p = 1) field distributions. Again, much of the fine storm structure 
is lost even in the highest frequency (5-mm band) channels. 

To reduce the complexity of the overall simulation, some potentially interesting aspects of 
the observation portion are tacitly avoided here. For instance, the polarization properties of the 
sensor and the upwelling brightness field are not considered. Also not considered are the effects 
of spot elongation or radiative transfer effects caused by different slant-path distances through 
the atmosphere and/or surface emisivity changes at off-nadir angles, as would be encountered 
at the extreme spots of LEO and wide-FOV GEO systems. These features are not considered 
to be important in the initial analysis of the expected performance of spatial resolution limited 
sensors. We note, however, that gain patterns for radiometric interferometers (or Synthetic 
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Aperture Interferometric Radiometers, SAIR’s [28]) can be analyzed with the same formalism 
and software as presented here for filled-aperture imaging radiometers. For the SAIR, the 
sensitivity pattern of the array would be substituted for the aperture gain pattern. Two rmnor 
differences are (1) the SAIR gain function is not positive definite, and (2) the influence of grating 
lobes cannot be neglected, especially in thinned-aperture SAIR arrays. 3 

Upon integration into the multispectral deconvolution experiments and precipitation re- 
trieval simulations, some issues that will be able to be quantitatively resolved using the above 
formalism include the tradeoff between beam efficiency and spot size, the effects of antenna gam 
pattern uncertainties and the effects of 3-dB raster sampling as opposed to Nyquist sampling. 


‘Due to the structural complexities oT^led- aperture sensor., SAIR*. may be the only practical method of 
implementing particularly large apertures. 
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4. Karhunen-Loeve Analysis 


The convolved Ta imagery has been analyzed for loss in observational degrees of freedom 
using the Karhunen-Loeve (KL) transformation (also called principal components analysis, or 
the method of empirical orthogonal eigenfunctions). The unitary KL transform rank-orders the 
dominant spectral modes in the T a imagery by the modes’ signal- to-noise ratios (SNR’s) (29, 30). 
The SNR’s of the Ta KL spectral modes are then compared to those of the Tb KL modes to 
determine the number of observational degrees of freedom in the observed imagery relative to 
the number of degrees in the full resolution imagery. 

The KL transform is determined by the eigenvectors of the spectral covariance matrix 
Ht a Ta • One way to estimate this matrix is to assume ergodicity, that is, that spatial averaging 
over the pixels of an image yields statistics identical to those obtained by ensemble averaging 
over pixels chosen from a large set of independent images. Strictly, this is only true for two- 
dimensionally ergodic random processes. In the simulated Tb and Ta imagery, this condition 
is not quite satisfied. However, the lack of a large ensemble of independent images requires this 
practice; the consequences of its effects will be discussed later in this section. 

Assuming spatial ergodicity, the spectral covariance for Ta becomes. 

%t a T a = (ft a - (Ta))(Ta ~ Ta))*) 


£*max 


sar it,,, 

‘maxjmax _ 1 


(H) 


where the outer expectation operator (•) is computed by space-averaging over the ensemble of 
all imaxjmax pixels in the image and we have arbitrarily defined Tg** to be the expected^value 
of Tb in the absence of hydrometeors. The eigenvectors and associated eigenvalues of Ht a t a 
are found by di&gonaliz&tion: 


*Rt a t a = E 


Aj 0 

0 A 2 



A/v 


(18) 


where the rows of the KL transformation matrix T are the eigenvectors of 1 It a t a - These 
eigenvectors provide a measure of the contribution from each channel of the brightness imagery 

to each KL image: 

lAii-KTAii-*?') ( 1Q ) 


An important property of the KL imagery is that all KL channels are uncorrelated in the 
sense of Eq. 17. Note this does not imply that they are statistically independent. Rather, 
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. . .. . - r j p . statistics Indeed, the KL transform (as based 

thev are independent only in their second order statistics, ino , - ... 

on Eq 17) produces complete statistical independence only m the case where Ta P . 

Gaussian^ rZ - process Since Ta doe, no, satisfy -hi, condition, the join, moments of F for 

orders higher than two are generally nonzero. 

For each KL channel, the associated eigenvalue A n is its total expected energy, deter 

mine the fraction of this energy containing geophysical information (^n, no.«Mh. «g»v^« 
1„ be compared with the observational noise energy in ,h. channel. Proved that the sen 
““vitLs or. are identicai over all sensor channels, this same noise energy «,« .PP«r ■» 

KL channel. This is a consequence of the unitary property of the transform matra 
cane the signal-to-noise ratio (SNR) for each KL channel is: 

SNR. = < 2 °> 

°Tn 

In practice, some sensor channels wiU be more sensitive than others, thus th^s will 
be different. To accommodate observational noise level variations among sensor channels ,t ,s 
necessary to prewhiten the brightness imagery prior to compulation of the covanance 


r 


Ai] ~ 


l/v ®57 

0 

0 1 /v^t 7 


1/v^TW J 

r i/v/^7 

0 

0 l/y/ffi 7 

l 

l/v^fw J 


i/v^t 7 o 


0 1 


Ai] 


( 21 ) 




*Tbv 


jr 


( 22 ) 


i/v^7 

where 7V is the noise-whitened antenna temperature imagery. The purpose of ' 

" tit. noise component of the signa. T A to he white. Since f is urn, ary ,h.s nots 
remains white under the KL transformation. That is, the noise energy in Fa « uncorrelated 
and distributed equally among all KL channels. This would not follow if the noise componen 
Z r ' l colored. Prewhifening by Eq. 22 also scales the transformed observe, .onai no,se 
energies in all KL channels to unity. Thus, the KL-channel SNR s are: 


SNR n = A„ - 1 


(23) 


, , , . i _ c 7? , Effectively, prewhitening allows comparison of the 

where the A n ’s are the eigenvalues of R t > a t > a • Mtectiveiy, P 6 . . . 

information contained in the brightness channels on a “channel-no.se eqmvalen, bas.s, 
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By rank ordering the SNR’s of the KL channels the relative importance of each combi- 
nation of sensor channels (also referred to as a “KL mode” or “KL image”) can be assessed. 

For example, for the full- resolution six-channel MIMR imagery, the most dominant KL image 
(i.e., with the largest SNR, 28.40 dB) shows the overall shape of the GCE 4 cell extremely 
well (Fig. 10a, leftmost image). The second and third modes (next two from left, 18.60 and 
4.80 dB, respectively) respond primarily to the bimodal brightness signature near the partially 
transparent edges of the cell. The other three modes are primarily noise (rightmost images), 
as suggested by the their low SNR’s. The adjoining eigenvector plots (Fig. 11a) reveal that 
the most dominant mode is essentially the spectral difference between the 23.8- and 89-GHz 
channels. We interpret this mode as responding to the strongly-varying scattering signature 
caused by ice aloft. The second KL mode is a combination of the scattering signature from the 
high frequency channels and the spectral difference between the 6- and 10.69-GHz channels, and 
appears to be sensitive to the radiometrically-thin periphery of the storm. The other modes are 
successively more difficult to interpret. Indeed, there is no guarantee that any KL mode has a 
meaningful geophysical interpretation. 

The KL-mode SNR’s for the full-resolution (7 b) and convolved (7 a ) brightness imagery 
can be used to estimate the loss in number of observable degrees of freedom caused by the broad- 
ness of the antenna pattern. As seen in Figs. 10b and c, the KL modes for the convolved MIMR 
imagery exhibit monotonically decreasing SNR’s, indicative of a progressive loss of spectral vari- 
ance caused by the spatial filtering of the antenna. A comparison of the SNR’s (Figs. 10 and 12) 
for this imagery shows a reduction of ~ 1.6 dB in the third mode from the full resolution to the 
4.4-m case, then additional reductions of ~ 1.3 dB and ~ 1.5 dB in the second and third modes 
(respectively) from the 4.4-m case to the 1.6-m case. These reductions are attributed to the 
relatively large dependence of modes 2 and 3 on the lower frequency channels for which spatial 
resolution is adversely affected by aperture size reductions. In contract, the dominant mode 
SNR’s show insignificant reductions with aperture size, reflecting the fact that the dominant KL 
mode depends mostly on the high frequency channels for which spatial resolution is better. 

Further insight into the effects of varying spatial resolution can be seen by examining 
' the eigenvectors in Fig. lla-c. Foremost, there is relatively little difference in the eigenvectors I 
for a given mode as the resolution is changed. This suggests that no gross changes in the 
spectral characteristics of the brightness imagery occur over the range of spatial resolutions 
being considered. However, some changes in the relative contributions from the six MIMR 
channels can be seen. In mode 1, the contribution from the 6- and 10.69-GHz channels decreases 
by a factor of ~ 2 from full- resolution to the 1.6-m case, thus indicating a reduction in signal 
variance in the low-frequency channels. However, the high frequency channels remain strong 
contributors, indicating that spatial resolution is good for all cases. In mode 2, the contribution 
from the low-frequency channels is relatively small, and changes sign from the full resolution to 
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the 1.6-m case. TO. suggests that the contribution to the bimodd structure of Kl mode 2 rom 
the lowest frequency channel is reduced by spatial filtering. In mode 3, the contribution from 
the low frequency channels decreases markedly from full resolution to 1.6-m resolution, again 

indicating a decreasing variance. 

However, caution must be exercised in interpreting the KL mode SNR’s in an absolute 
„ The first (and most obvious) reason is that the SNR’s have been derived from only a 
tingle time frame (GCE 4) of a computer-generated brightness held. Ideally, a large ensem e 
of real Image, would be used. A more subtle reason arise, from the non-ergodicity of the 
precipitation cell brightness temperatures. Since each cell is Unite in extent, an essentia^ question 
Li«. as to what fraction of the image area should be padded with noisy clear-ai, bngh.n«s 
temperatures. The larger this padding, the smaller the multing eigenvalues V Thus he 
eigenvalues (and hence the SNR’s) depend on the size of the cell in relation to the size o e 
image. To quantify this relation, let the number of pixels used in the space-averaging process 
of Eq. 17 be A, = imaxjmax. end the resulting eigenvalues mid SNR’s >>«*.“'» S &>• 
respectively. Then padding the image with additional noisy, dear-air pixels will yield a larger 
image with A" = -Wlmax Pi«ls, where N' r > A’,. The eigenvd.es of the padded imagery 
will be correspondingly smaller, as determined by the number of additional pixels. 

v _ , *£ + (24) 

A„ - a„ A ,, + w . 


Then from Eq. 23, it follows that: 


SNR' = SNR 




( 25 ) 


or, in decibels: /n\ 

SNR^(dB) = SNRn (dB) + 101og 10 I (26) 

Thus, the KL mode SNR’s for precipitation cell imagery are mearungful only to within an 
arbitrary additive constant which depends on the size of the dew-air portion outside the cell 
However, for identically sized images, the SNR’s can be meaningfully intercompared. Visual 
inspection of the KL imagery can dm provide a meaningful quditative determination of the 

usefulness of a mode. 

For the nine-channel GEO LSA system there are nine KL modes, of which perhaps 4-6 
provide significant information (Figs. 13a-c). Agdn. the KL mode, for the convolved imagery 
exhibit monotonicdly decreasing SNR’s, indicating a progressive loss of spectrd variance as the 
aperture size decreases. A comparison of the SNR’s (Fig. 15) show, a reduction of ~ 2 in 
the dominant mode SNR from the full resolution to either the 15-m or 40-m cases, and much 
larger reductions (~ 10 dB or more) for the second and higher order modes. 
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Again, the eigenvectors for a given mode (Figs. 14a-c) exhibit no gross changes as the 
spatial resolution is varied. As in the MIMR case, mode 1 consists primarily of the difference 
between the high frequency channels (which cool in response to scattering by ice aloft) and 
the low frequency channels (which warm in response to rain and cloud water). This effectively 
produces a map of the cell. Thus, the mode 1 eigenvector is a good indicator of the relative 
importance of the high-frequency millimeter and submillimeter channels (e.g., 89, 166, 220, 340 
and 410 GHz) for precipitation cell mapping. While the variances of the channels above 89 GHz 
fall off with frequency, the diffraction-limited resolution of these channels increases. From the 
KL analysis, it appears that this tradeoff results in the the high frequency channels contributing 
most of the variance to the dominant mode. As aperture size (and hence resolution) decreases, 
the contribution from the low frequency channels further diminishes, illustrating the loss of 
variance from these channels caused by spatial filtering. 

For the eight channel GEO LSA system, approximately three KL modes provide signifi- 
cant information (Figs. 16a-c). As seen in the eigenvectors for these channels (Figs. 17a-c), the 
5-mm band oxygen channels respond in approximately the same manner to precipitation. Of 
course, the primary purpose for the five oxygen channels would be for temperature sounding. If 
precipitation were the only concern, only the 50.3 GHz channel would be used. The dominant 
KL mode is primarily comprised of the difference between those channels that respond to pre- 
cipitation by primarily increasing their brightness and those that respond by decreasing their 
brightness. As the aperture size decreases, the SNR of the dominant mode drops by 2 dB for the 
40-m case, then another 2 dB for the 15-m case(Fig. 18). The SNR for all higher order modes 
are more significantly affected. 

It is of interest to note that the eigenvectors for a given system (as computed using Eq. 18) 
are relatively insensitive to the time of sampling, that is, all frames catalogued in Table 1 result in 
similar sets of eigenvectors. (The SNR’s differ, but this is the result of the size and intensity of the 
precipitating region in relation to the clear air surrounding it.) Thus, the eigenvectors for a given 
channel set can be considered to be quite “universal” in that they are unique to the system as 
well as to the class of meteorological event (in this case, convective precipitation). This has also 
been observed in 118-GHz imagery of storm cells, where it has been found that approximately 
2-3 unique observable modes are available from realistic 118-GHz imaging systems observing 
convective precipitation (23). This suggests that retrieval algorithms can be constructed using a 
reduced set of inputs variables, namely, the KL modes. In addition, it suggests that compression 
of images for coding and archival purposes can be accomplished with a high ratio of compression 
(typically a factor of five to ten). 
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5. Multispectra] Antenna Pattern Deconvolution 


The KL analysis provides a measure of the number of spectral degrees of freedom pro- 
duced by precipitation in the observed T A imagery. It does not, however, provide a measure of 
the spatial-frequency information lost by filtering the Tb imagery with the multispectral sensi- 
tivity pattern. To quantify the loss of spatial-frequency information and at, the same time, to 
demonstrate the usefulness of the KL spectral modes, a statistically optimal deconvolution of the 
multispectral spot pattern and the associated error caused in the observation and deconvolution 
process is being investigated. The multispectral statistical deconvolution (MSD) operator is 
necessarily linear, and is formulated to minimize the expected error between the full-resolution 
brightness imagery and the deconvolved imagery. 

The basic procedure for MSD, described by Rosenkranz (31^ requires convolving the 
observed T, t with an optimal shift-invariant deconvolution operator 3 ,j to obtain an estimate 

Tb of Tb : ^ 

Tb < ; = 

By requiring that the error between Tb and Tb be uncorrelated with the observations T A . 

{(T B ij - Tb^T'a,?) = 0 ( 28 ) 

where all possible values of « are considered, the minimum mean-square-error estimator 
5 can be shown to follow: 

i,i ••T, - <T,« -TU-j) <»> 

This is recognized as an application of the orthogonality principle in estimation theory. In con- 
trast to the spatial averaging used in the KL analysis, the ensemble averaging must be performed 
over a number of statistically independent images. The reason for the difference is that spatial 
correlations (i.e., correlations between different pixels) are required to be incorporated into the 
statistics in Eq. 28, whereas the KL analysis required no inter-pixel statistical information. 

Using the two-dimensional discrete Fourier transform (denoted by T), 3 can be repre- 
sented in the discrete spatial-frequency domain as: 

, •max -1 imax -1 _ • , ... , 

^ _ l y y e -j2*(««/‘max + vj/Jmax) (30) 

*maxJmax , m0 Ja0 

where the integer variables u,t> denote discrete spatial frequencies separated by increments of 
(imaxAJl)- 1 and (jmaxAfl) -1 , respectively. Using the convolution-multiplication property of 
Fourier transform pairs and Eq. 14, we can write D as: 

s „„ = = ^[( t b „ « n .,,.,) i ( ri ( T <0 •. tv ,.,)])- 1 
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= F\{ Tb i) **{Ui) **T B -i,-i + W-i.-i )* ) 1 

= ?[(Tbh **($,} **T b -,,->)')] 

(/■[( **T B ij ) **T B -i,-j y ) + ( *0 "rf-i.-j )]) 

= ) t$„\ ( m,\' < TV*-,) Ws.,)' > m,i + (m/i^]’)r’ 

= 1>TbTb (S tiV 2 'TbTb uv^v + 5im#v) 


(31) 


where is the transform of the gain matrix f,j, 7>t b T b »v * 8 the^multidimensional power 
spectral density matrix for the brightness temperature vector Tb, J»n«* is the multidimen- 
sional power spectral density matrix for the observation noise vector IT, and * is the Hermitian 
transpose. Since IT is white, we have: 


^nniiv “ i;] “ 


1 


•maxJmax 


4, 0 

0 


2 

a T7 


a \ N J 


(32) 


which is identical for all spatial frequencies u, v. The second-to-last step in Eq. 31 uses the sta- 
tistical independence of the brightness process T B and observation noise Tl. Note that *§t b t b uv 
contains all information on both the spatial and intrachannel spectral statistics of the random 
brightness process. Here, the use of the multidimensional power spectral density matrix implic- 
itly assumes that the brightness random process is wide-sense stationary. This is a slightly less 
restrictive assumption than the ergodic assumption used in the KL analysis. 


Thus, the optimal linear deconvolution becomes: 

fauv = V uv rF Aij ] (33) 

with 

Tea = r-'PB-) ( 34 ) 

where ^""Ms the two-dimensional discrete inverse Fourier transform. This is equivalent to Eq. 27, 
except written in the spectral domain. Defining the MSD error to be: 

= Tb.j - Tbo (33) 


The associated MSD spectral error covariance is [31]: 

= ^T B T B uv ~ ^T B T B uv^uv((j U v^T B T B uv^w + ^nnuv) l ^lv^T B T B uv 

(36) 

Note that this error is always less than the a-priori power spectral density for Tb, illustrating 
that the optimal deconvolution operator can only reduce the error spectral density. The spectral 
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domain MSD operator uses spatial information from the highest resolution channels to improve 
the deconvolution of the lower resolution channels. Such improvements are possible to the 
extent that the brightnesses between the low frequency and high frequency channels exhibit 
known correlations [32, 33, 34]. These correlations can be due to either horizontal rain cell 
structure or spectral scattering or absorption similarities. 


The effects of antenna gain uncertainties are manifested as additional non-white obser- 
vation noise. To see this, consider the gain matrix to be randomly perturbed by a random gain 


error 


Ta: 



(37) 


where is the assumed gain matrix from (for example, from Eq. 7). The observed antenna 
temperature imagery now becomes: 


T AiJ = **T b , } + ? **Tbh + 

= V-i'-i **Tbh + 


(38) 


where It' includes contributions from both radiometric integration noise and gain pattern un- 
certainties. The noise covariance matrix for if becomes: 

En'n'.j = " >') = ^ Bi ' )) 

(39) 


<yfj 0 
0 a\ 2 


°TN J 




(40) 


where the integration noise is assumed to be uncorrelated with the gain uncertainty and the 
averaging in the second line of Eq. 40 is taken over an ensemble of antenna gain perturbations. 
In the spectral domain, the corresponding noise power spectral density matrix is: 

7n'«' *v = 2„„ $ T B T B uvttuv) + ( ^00 ^°° ) 6 « 

= 

= 5 »"” + £ r .„ + (41) 
where the matrix f£ T , w is Hermitian square root of the brightness power spectral density 

7*7® uv ^ =*1/2 1 ^1/2 (42) 

■->7*7* uv ~ ^ TbTb vx s T b Tb*v v ' 
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The first term in Eq. 41 is the flat receiver noise power density spectrum (Eq. 32). The second 
term describes additional noise caused by brightness fluctuations due to clouds and raincells as 
well as surface features such as shorelines, small islands and ice floes. This portion of the noise 
spectrum is a low pass process, but is assumed to be zero mean so that it provides virtually no 
contribution at zero spatial frequency (u = 0,t> = 0). The third term describes uncertainty in 
the offsets of the multispectral imagery caused by uncalibrated spillover of the antenna beam 
into cold space. This is a zero-frequency or DC contribution to the noise power spectrum. 

Presumably, most gain perturbations would be caused by deformations in the reflector 
and/or supporting structure, and hence would affect all channels in a correlated manner. More- 
over, the gain perturbations within any given channel are not expected to be independent over 
all pixel values »' and j. Rather, whole sidelobes encompassing many pixels might increase and 
decrease simultaneously. Thus, the matrices developed in Eqs. 40 and 41 cannot be considered 
to be diagonal, nor can it be considered to be uncorrelated among pixels. 

Using Eqs. 36 and 41, the effects of antenna gain errors caused by both structural pertur- 
bations and uncalibrated spillover can be analyzed for their effects on the optimum deconvolved 
brightness maps. However, to examine their effects on the retrieved precipitation parameters, 
the end-to-end retrieval problem simulation outlined in the next section must be considered. 

The MSD operator can also be developed as an optimal linear estimator in the spatial 

domain: __ 

f B =Tf A (43) 

where 7 a a vector consisting of all pixels from all channels having any significant correlations 
with the brightness at the location of interest. The determination matrix if is given by: 

T? - Ht b T a (Bt a T a +^fin) -1 (44) 

H a 0 is the joint spatial-spectral covariance matrix between the multispectral vector signals o and 
/?. It is noted, however, that this method results in particularly large determination matrices. 
An advantage is that the brightness imagery need not considered to be stationary, that is, the 
statistics can be allowed to vary from point to point within the image. This would be useful for 
analysis of imagery observed near boundaries such as coastlines, fronts or image edges. 

Currently, demonstrations of the spectral and spatial domain MSD operators using the 
GCE data are being constructed. Of particular interest will be the effect of the noise-resolution 
tradeoff first described by Backus and Gilbert (35). Indeed, as the resolution AR of the decon- 
volved map is reduced, the noise at each deconvolved pixel increases. Note that this tradeoff will 
be inconsequential for area-averaged linear parameter retrievals since the noise is averaged out 
over the image again in the end. However, for non-linear retrievals, there may be some optimal 
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deconvolution grid spacing for a given system. The high resolution brightness maps and the 
intermediate MSD step will allow this search to be performed. 
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(i. Precipitation Retrieval Simulation 


This component of the PM MRS will demonstrate simulated precipitation retrieval algo- 
rithms for multispectral low-pass filtered microwave imagery. Initial retrievals of precipitation 
will use the single-step nonlinear statistical technique demonstrated in (23). This technique is 
based on a non-linear mapping of the Karhunen-Loeve image data Tc A ij to various meteorological 
parameters, for example, surface rain rate, integrated water content, integrated ice content and 
mean ice particle size. The non-linear mapping is used to estimate a value for each of these 

parameters at each pixel. 

It is expected that more accurate retrieval methods will employ a statistical iterative 
technique. This scheme has proved to be useful in retrieving relative humidity profiles [36, 37]. 
A statistical iterative retrieval technique based on the incrementally linear relationship between 
small changes in brightness temperature and small changes in surface precipitation rate (via 
the forward raditive transfer process) is planned for implementation in the PMMRS [38]. This 
technique is similar to the linear spectral domain MSD technique except that ~D UV becomes an 
incremental predictor matrix, t>3 uv , and STbuv and 6R UV (where 6R = T\br) is the incremental 
change in the rain rate field r„) replace T Auv and t fluw , respectively, in an iterated Unear 
statistical operator. This relationship, similar to (Eq. 43) is given by 

6R UV = 6T A uv ( 45 ) 

where fT Auv is the error between the observed antenna temperature and the computed antenna 
temperature: _ 

6T Auv = T Avv - ] ( 4fi ) 

In the above, we assume that Tsij is a known (nonlinear) function of r,j. A requirement is that 
a quick and accurate method of computing the forward transfer problem be available. This may 
be based on a piecewise-linear approximation to the forward radiative transfer relationship. 

To take full advantage of the spatial resolution available from high-frequency channels, it 
may seem that multispectral deconvolution of the antenna gain pattern must be first performed. 
However, the deconvolution aspect of the problem can be integrated into the nonlinear precip- 
itation parameter retrieval using a spectral-domain nonlinear estimator. In these schemes, the 
effects of antenna pattern errors will be examined by their effect on the modified instrument 
noise covariance (Eq. 41). Iteration ceases when the &T Auv is zero to within the noise of the 
sensor, as determined by l n >n>- To reduce the complexity of the calculation, the KL modes can 

be used. 
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7. Discussion 


At this time, the microphysical cloud and raincell data compilation, forward radiative 
transfer calculations, sensor observation simulations, and Karhunen-Loeve transformation and 
degree of freedom calculations have been implemented. Although we cannot yet quantify the 
available retrieval accuracy of candidate systems, the loss in number of observable degrees of 
freedom as aperture size is reduced is readily seen. The effect of this information loss will be more 
fully understood upon completion of the end-to-end retrieval simulation. Future implementation 
of the deconvolution and parameter retrieval components of the PMMRS will improve the ability 
to determine (quantitatively) the available accuracy of existing and future spaceborne passive 
microwave meteorological systems and to design optimized sensor-specific multispectral retrieval 

algorithms. 
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Tkble 1: Catalogue of synthesized brightness imagery used in the numerical simulations. 
“ H Horizontal Image 

Brightness # Resolution Size 


Map Set 

Channels 

A R (km) 

(km) 

Comments 

GCE 1 
GCE 2 
GCE 3 
GCE 4 
GCE 5 
CP-2 

■ 

m 

96x96 

96x96 

96x96 

96x96 

96x96 

41x41 

Tao & Simpson GCE data, t — 126 min 
Tao & Simpson GCE data, t = 138 min 
Tao & Simpson GCE data, t = 174 min 
Tao & Simpson GCE data, t = 210 min 
Tao & Simpson GCE data, t = 234 min 
COHMEX CP-2 data, July 11, 1986 


(window channels only) 


Table 2: Scaling factor a (p) and a perture efficiencies for v arious aperture taper parameters p. 

T Aperture 
Efficiency 
p a r) A (%) 

0 I 1.02 I 100 

1 1.27 75 

2 1.47 56 

3 1.65 44 

4 1.81 36 
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Table 3: System sensitivity parameters for the MIMR simulations, 
based on the single-spot integration time for sampling at the Nyquist 

frequency channel. — — — . r - 

Channel Frequency Bandwidth 7*y, Integration 
# (GHz) (MHz) 1 (K) 1 Time (msec) | 


The integration time is 
resolution of the highest 


= 1.89 


6.00 

20 

250 

10.69 

20 

350 

18.70 

200 

400 

23.80 

400 

400 

36.50 

1000 

500 

89.00 

Kmili 

650 


0.030 0.32 

0.16 

0.45 

0.22 

0.16 

0.079 

0.12 

0.059 

0.10 

0.049 

0.05 

0.024 


?s * 0.27 km, %dB.N 
I 250 I 0.0039 I 


= 0.69 

0l8~' 

1.23 

0.44 

0.33 

0.27 

0.14 


0.16 

0.22 

0.079 

0.059 

0.049 

0.024 
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Table 4: System sensitivity parameters for the LSA simulations. The integration time is based on 
the single-spot integration time for sampling at the Nyquist resolution of the highest frequenc) 


channel. 


Channel 

# 

Frequency 

(GHz) 

Bandwidth 

(MHz) 

T. v . 

(K) 

Integration 
Time (msec) 

&SS n 

(K) 

OTn 

(K) 

LSA: nine channels, D 

= 15 m, R s = 

= 0.88 km, %dB,N = 

2.22 ki 

m 

1 

6.00 

20 


250 

1.2 

1.63 

0.95 

2 

10.69 

20 


350 


2.29 

1.33 

3 [ 

18.70 

200 


400 


0.83 

0.48 

4 

36.50 

1000 


500 


0.46 

0.27 

5 

89.00 

6000 


650 


0.24 

0.14 

6 

166.00 

4000 


1000 


0.46 

0.27 

7 

220.00 

3000 


2000 


1.07 

0.62 

8 

340.00 i 

3000 


3000 


1.60 

0.93 

9 

410.00 

3000 


5000 


2.67 

1.56 

LSA: nine channels, D 

= 40 m, Rs 

= 0.33 km, %d b,N = 

= 0.83 km 

1 

6.00 

20 


250 

0.16 

4.34 

0.95 

2 

10.69 

20 


350 


6.10 

1.33 

3 

18.70 

200 


400 


2.21 

0.48 

4 

36.50 

1000 


500 


1.22 

0.27 

5 

89.00 

6000 


650 


0.64 

0.14 

6 

166.00 

4000 


1000 


1.22 

0.27 

7 

220.00 

3000 


2000 


2.85 

0.62 

8 

340.00 

3000 


3000 


4.26 

0.93 

9 

410.00 

3000 


5000 


7.11 

1.56 

LSA: eigl 

it channels, D — 15 m, 

Rs 

= 6.53 km, i?3-dB,N 

= 16.6 

km 

1 

18.700 

200 


400 

65.3 

0.11 

0.48 

2 

23.800 

400 


400 


0.078 

0.34 

3 

36.500 

1000 


500 


0.062 

0.27 

4 

50.300 

180 


550 


0.16 

0.70 

5 

52.800 

400 


550 


0.11 

0.48 

6 

53.596 

170 


550 


0.16 

0.70 

7 

54.400 

400 


550 


0.11 

0.48 

8 

54.940 

400 


550 


0.11 

0.48 

LSA: eight channels, D = 40 m, 

Rs 

= 2.45 km, %dB,N 

SJ 

to 

II 

km 

1 

18.700 

200 


400 

9.18 

0.29 

0.48 

2 

23.800 

400 


400 


0.21 

0.34 

3 

36.500 

1000 


500 


0.16 

0.27 

4 

50.300 

180 


550 


0.43 

0.70 

5 

52.800 

400 


550 


0.29 

0.48 

6 

53.596 

170 


550 


0.43 

0.70 

7 

54.400 

400 


550 


0.29 

0.48 

8 

54.940 

400 


550 


0.29 

0.48 
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Figure 1: Schematic diagram of the PMMRS. 


Figure 2: Computed window-channel brightness maps for three GCE simulations: (a) GCE 2, 
(b) GCE 3, (c) GCE 4, and for the COHMEX CP-2 data (d). The channel frequencies (in GHz) 

are indicated. 

Figure 3: Computed brightness maps for channels near the microwave oxygen lines for GCE 3. 
The channel frequencies (in GHz) are indicated. 

Figure 4: Computed brightness maps for channels near the microwave water vapor lines for 
GCE 3. The channel frequencies (in GHz) are indicated. 


Figure 5: Scan patterns of canonical (a) geosynchronous and (b) low-Earth orbit imaging sys- 
tems. 


Figure 6: Energy distribution functions for linearly-polarized circular apertures of varying am 
plitude tapers p, indicating efficiencies at the 3-dB beamwidth (□) and first null (•). 


Figure 7: M1MR brightness temperature imagery, as would be observed from low Earth orbit: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 4.4-m 
quadratically-tapered aperture distribution, respectively. The 3-dB spot sizes (in km, FWH ) 
are indicated for each frequency in (b) and (c). 


Figure 8: Representative LSA brightness temperature imagery for nine window channels, as 
would be observed from geosynchronous orbit: (a) full resolution, and (b,c) after convolution 
with the gain patterns of circular 40-m and 15-m quadratically-tapered aperture field distribu- 
tions, respectively. The 3-dB spot sizes (in km, FWHM) are indicated for each frequency in (b) 

and (c). 

Figure 9: Representative LSA brightness temperature imagery for eight low-frequency channels, 
as would be observed from geosynchronous orbit: (a) full resolution, and (b,c) after convolution 
with the gain patterns of circular 40-m and 15-m quadratically-tapered aperture field distribu- 
tions, respectively. The 3-dB spot sizes (in km, FWHM) are indicated for each frequency in (b) 

and (c). 

Figure 10: Rank-ordered KL mode imagery for the six-channel MIMR system: (a) full resolution, 
(b,c) after convolution with the gain pattern of a circular 4.4-m and 1.6-m quadratically-tapered 
aperture distribution, respectively. The five most dominant modes are shown with the respective 
SNR’s indicated above each image. 

Figure 11: Three most dominant rank-ordered eigenvectors for the six-channel MIMR system: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 4.4-m and 1.6-m 
quadratically-tapered aperture distribution, respectively. 

Figure 12: SNR ranking for the full-resolution and convolved MIMR imagery. 
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Figure 13: Rank-ordered KL mode imagery for the LSA nine-channel GEO system: (a) full res- 
olution, (b,c) after convolution with the gain pattern of a circular 40-m and 15-m quadratically- 
tapered aperture distribution, respectively. The five most dominant modes are shown with the 
respective SNR's indicated above each image. 


Figure 14: Three most dominant rank-ordered eigenvectors for the nine-channel LSA system: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 40-m and 15-m 
quadratically-tapered aperture distribution, respectively. 


Figure 15: SNR ranking for the full-resolution and convolved nine-channel LSA imagery. 


Figure 16: Rank-ordered KL mode imagery for the LSA eight-channel GEO system: (a) full res- 
olution, (b,c) after convolution with the gain pattern of a circular 40-m and 15-m quadratically- 
tapered aperture distribution, respectively. The five most dominant modes are shown with the 
respective SNR’s indicated above each image. 


Figure 17: Three most dominant rank-ordered eigenvectors for the eight-channel LSA system: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 40-m and 15-m 
quadratically-tapered aperture distribution, respectively. 


Figure 18: SNR ranking for the full- resolution and convolved eight-channel LSA imagery. 
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Figure 1: Schematic diagram of the PMMRS. 
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Figure 2: Computed window-channel brightness maps for three GCE simulations: (a) GCE 2, 
(h) C.CE 3, (c) GCE 4, and for the COIIMEX CP 2 data (d). The channel frequencies (in C.llz) 
are indicated. 
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Figure 2 (H) 































Beam Efficiency Vb(9) 


Beam Efficiency vs. Angle 
E=E 0 (1-(2p/D) 2 ) p with D=30A 



Beam Angle Parameter U=(7TD/A)sin0 


Figure 6: Energy distribution functions for linearly-polarized circular apertures of varying am 
plitude tapers p, indicating efficiencies at the 3-dB beamwidth (■) and first null (•). 
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Figure 7: MIMR brightness temperature imagery, as would be observed from low Earth orbit: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 4.4-m and 1.6-m 
quadratically-tapered aperture distribution, respectively. The 3-dB spot sizes (in km, FWIIM) 
arc indicated for each frequency in (b) and (c). 
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Figure 8 (b) 
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Figure 9: Representative I,SA brightness temperature imagery for eight low-frequency channels, 
as would be observed from geosynchronous orbit: (a) full resolution, and (b,c) after convolution 
with the gain patterns of circular 40-m and 15-ni quadratically-tapered aperture field distribu- 
tions, respectively. The 3-dH spot sizes (in km, FWIIM) are indicated for each frequency in (b) 















Figure 9 (c) 
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Figure 11: Three most dominant rank-ordered eigenvectors for the six-channel MIMR system: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 4.4-m and 1.6-m 
quadratically-tapered aperture distribution, respectively. 
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Figure 12: SNR ranking for the full-resolution and convolved MIMR imagery. 
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Figure 13: Rank-ordered KL mode imagery for the LSA nine-channel GEO system: (a) full res- 
olution (b,c) after convolution with the gain pattern of a circular 40-m and 15-m quadratically- 
taperrd aperture distribution, respectively. The five most dominant modes are shown with the 
respective SNR *s indicated above each image. 
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Figure 15: SNR ranking for the full-resolution and convolved nine-channel LSA imagery. 
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Figure 17: Three most dominant rank-ordered eigenvectors for the eight-channel LSA system: 
(a) full resolution, (b,c) after convolution with the gain pattern of a circular 40-m and 15-m 
quadratically- tapered aperture distribution, respectively. 
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Figure 18: SNR ranking for the full-resolution and convolved eight-channel LSA imagery. 
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